Gravitational waves from hyperbolic encounters 
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The emission of gravitational waves from a system of massive objects interacting on hyperbolic 
orbits is studied in the quadrupole approximation. Analytic expressions are derived for the gravita- 
tional radiation luminosity, the total energy output and the gravitational radiation amplitude. An 
estimation of the expected number of events towards different targets (i.e. globular clusters and the 
center of the Galaxy) is also given. In particular, for a dense stellar cluster at the galactic center, a 
| rate up to one event per year is obtained. 
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General Relativity predicts that a system of interacting massive objects emits gravitational waves which propagate 
' in the vacuum with the speed of light. A lot of studies have been devoted to the description of gravitational radiation 
emitted by systems of two interacting stars where amplitude, power and all physical quantities of such a radiation 
i— - n ' strictly depend on the configuration and the dynamics of the system. The seminal papers by Peters and Mathews 

0. 01 , which investigate the gravitational wave emission by a binary system of stars (on circular or elliptic orbits) in 
the quadrupole approximation, have been extended in several directions (see e.g. Q and references therein) in which 

}— i the problem is studied by both analytical and numerical approaches. 

On the other hand, depending on their approaching energy, stars may interact also on unbound orbits (parabolic or 
hyperbolic) and, in this case, one expects that gravitational waves are emitted with a peculiar signature of a "burst" 
£T) i wave-form having a maximum in correspondence of the peri-astron distance. The problem of Bremsstrahlung-like 
gravitational wave emission has been studied in detail by Kovacs and Thorne [H by considering stars interacting on 
unbound orbits. 

Here, we face this problem discussing the dynamics of such a phenomenon in the simplest case where stars interact 
on unbound orbits and the peri-astron distance is much larger than the Schwarzschild radius of the stars. In this 
case, the quadrupole approximation holds and useful analytic quantities (as the energy emitted by the system per- 
unit time and gravitational wave amplitude) can be derived. 

In this letter, after reviewing the main features of hyperbolic encounters between two stars, the emission of grav- 
itational waves is studied in the framework of the quadrupole approximation. Particular attention is devoted to the 
gravitational radiation luminosity. Then, we derive the expected gravitational wave-form and discuss the detection of 
such events which could greatly improve the statistics of possible gravitational wave sources. 

The study of gravitational wave emission from massive objects interacting on hyperbolic orbits can be started by 
analyzing the geometry of hyperbolic encounters. Without loss of generality, let us consider a mass Mi moving in the 
gravitational potential $ generated by a second mass Mi at rest in the center O of the reference frame (see Fig. [TJ. 
Let ON be a reference direction in the orbital plane. Then, the position of Mi (in P) is specified by the vector radius 
r and by the polar angle <p which the vector radius forms with ON, <j> being measured in the direction of the star 
motion. Obviously, both the vector radius and the polar angle depend on time as a consequence of the star motion, 

1. e. r = r(t) and <j> — 4>(t). With this choice, the velocity v of the mass Mi can be parameterized as 

v = v r r + v^cf) , (1) 
where the radial and the tangent components of the velocity are, respectively, 

dr d<j) 

Vr = Tt v * = r lTf (2) 

In this case, the total energy and the angular momentum, per unit mass, read out 
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Figure 1: The geometry of the hyperbolic encounter. The mass Mi (in P) is moving on hyperbolic orbit (continuous line) with 
focus in O where the mass M2 lies. The motion of Mi is described by the vector radius r and the polar angle 0. The vector 
radius ro (corresponding to the polar angle <f>o) represents the peri-astron distance, i.e. the distance of the closest approach 
between the two interacting stars. 



and 

respectively. At this point, it is useful to adopt the variable u = 1/r so that eq. can be rewritten as 

du\ 2 2 2$(r) 2E 

-j + . 2 + ^ = ^=const., (5) 
which, differentiated with respect to u, gives 

d 2 u 1 d$(u) 



L 2 du 



(6) 



The force acting on the unit mass is F(r) = = -GM 2 u 2 which can be expressed as the gradient of the 

gravitational potential <&. Hence, the last equation can be recast in the form [j| 

d2u 1 T ,f 1 \ 

and thus 



GM-2 
T 2 



whose general solution, giving the parameterized trajectory followed by the unit mass particle, is 

u($=tfoos(0-0o) + -£r • (9) 

Here, C is a constant depending on the initial conditions and 4>o is the polar angle corresponding to the peri-astron 
distance, i.e. the distance of closest approach between the two interacting particles. Note that after the interaction, 
the particle is deflected by the angle <fid — 2(f>Q. 
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For a system of two objects (of mass Mi and M 2 ) interacting on hyperbolic orbits, the two-body problem can be 
reduced to the problem of a reduced mass particle /1 = M 1 M 2 /(M 1 + M 2 ) moving in the gravitational field generated 
by a total fictitious mass M = Mi + M 2 . Let us assume that M 2 is at rest while Mi is moving with initial (at infinite 
distance) velocity "Uo with impact parameter b (see Fig. [1}. 

Considering what previously stated, it is then clear that the trajectory followed by the reduced mass particle is 

(a\ n (a j \ 1 G(Mi+M 2 ) 
u{4>) = C cos(0 - <f>a) H : — — '- . (10) 

Differentiating the previous relation with respect to the time t and turning back to the variable r, one obtains 

(IT 

_ = C*Lsin(0-0 o ) , (11) 
by which it possible to determine the value of the constant C through the initial conditions of the motion, i.e. 



which can be written, being L = bvo, as 



Lsm(4> ) 



sin <po 



Alternatively, the initial conditions of the motion can be used directly in eq. (TIT)]) . In this case, one finds 

b z cos <po 

which compared with eq. (|13[) gives the result 

bv 2 

ta ^-G(M^' (15) 

Clearly, eci. tfTB")) is meaningless for tfio — 90°, which means that, for this critical value, there is no interaction between 
the stars which are at an infinite distance each other. In other words, stars are very far when (f> approaches to 90°. 
Finally, the equation of the orbit followed by the reduced mass particle turns out to be Q 

b sin 4>n 

(16) 



COs(</> — (f>o) — COS(/>o ' 

which allows to determine the modulus of the radius vector r as a function of the polar angle <f> once the initial 
conditions are known. With this set of equations at hand, we can estimate both the gravitational wave luminosity 
of the system and the gravitational radiation wave-form in the quadrupole approximation. We are adopting such an 
approximation since we are considering "binary" systems whose distance is larger than the "capture" distance. In 
other words, we are considering situations and initial conditions where pointlike approximation of stars holds and the 
system remains unbounded. In such cases, quadrupole approximation works and results are reasonable 0]. 

The Einstein field equations give a description of how the curvature of space-time, at any event, is related to the 
energy-momentum distribution at that event. In the weak field approximation, it is found that systems of massive 
moving objects produce gravitational waves which propagate in the vacuum with the speed of light. It can be shown 
that the energy emitted per unit time, in the form of gravitational radiation (after integrating on all the gravitational 
wave polarization states), is [8J 



dE G(D$Dm 



dt 45c 5 



(17) 



where the dot represents the differentiation with respect to time, the symbol () indicates the scalar product and 
the quadrupole mass tensor is defined as 

A 3 =^m a (3xX-V!), (18) 
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r a being the modulus of the vector radius of the a — th particle. 

It is then possible to estimate the amount of energy emitted in the form of gravitational waves from a system of 
massive objects interacting on hyperbolic orbits. In this case, the components of the quadrupole mass tensor are 

D n = ^r 2 (3cos 2 0- 1) , 
L> 22 =M?~ 2 (3sm 2 0- 1) , 

(19) 

D12 = D 2 \ = 'ifir 2 cos 0sin , 
D33 = —fir 2 , 

which can be differentiated with respect to time as required in eq. (|17p . In doing this, we can use some useful relations 
derived above, in particular eqs. ([3]), (jT5]) and (fTB|) . It is straightforward to show that 

D W D <M = MW miM , (20) 



where /(0, fa) is given by 



/(0, fa) = sin 4 (0 o - 0/2) sin 4 (0/2) tan" 2 O sin" 6 O 

x [150 + 72 cos 20o + 66 cos 2(0o-0) (21) 
-144 cos (20o - 0) - 144 cos 0] . 



Hence, the energy emitted by the system per unit time is 

dE _ Z2GL Q n 2 
~dt ~ 45c 5 6 8 



/(0,0o), (22) 



or, equivalently, 



which, for Mi = M 2 , can be rewritten as 



dE 32G^ V 2 

dF = -^5^ / ^^ o) ' (23) 



^ = -W /( ^ 0) ' (24) 

r s being the Schwarzschild radius of the mass. 

The total energy emitted in the form of gravitational radiation during the interaction is given by 

f°° rip 1 

AE= / \ — \dt . (25) 



o 



dt 



Since eq. ^ holds, we can adopt the angle as a convenient integration variable. In this case, the energy emitted 
for 0! < < 02 is 

a 771/ j. j \ ±r s mvl sin 2 0o /(0,0 O ) 

A-E(0i,0 2 ) = „ / jj — t \ -TTn d(f> , (26) 

45c J & J fa [cos (0 - 0o ) - cos O ] 2 

and the total energy can be determined from the previous relation in the limits fa — > and fa — > 20o. Thus, one has 

A£=%^F(Mo), (27) 
where F(b,vo) only depends on the initial conditions and it is given by 

F(b,v ) = [7206 tan 2 O sin 4 o ]~ 1 x (26280 o 

+232800 cos 20o + 1440 o cos 40 o (28) 
-1948 sin20 o - 301 sin 40 o ) . 
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In other words, the gravitational radiation luminosity strictly depends on the configuration and kinematics of the 
binary system and it improves at short b and high vq. 

Direct signatures of gravitational radiation are its amplitude and its wave-form. In other words, the identification of 
a gravitational radiation signal is strictly related to the accurate selection of the shape of wave- forms by interferometers 
or any possible detection tool. Such an achievement could give information on the nature of the gravitational wave 
source, on the propagating medium, and even, in general, on the gravitational theory producing such a radiation 
9, 10]. It is well known that the amplitude of the gravitational waves can be evaluated as 

Of 

h jk (t,R) = j^& k , (29) 

R being the distance between the source and the observer and {j, k] = 1, 2. 

Considering our binary system and the single components of eq. (|29| , it is straightforward to show that 

h u = as c S c 2 (0 o )M^ [ 13 cos _ 12 cos 2< f, + 3 C os 30 

-2 cos(0 - 20 o ) + 3 cos(30 - 20 o ) - 12 cos(20 o ) 

-6 cos 2(0o - 0) - 6 cos 2(0 + <j> ) + 15cos(0 + 20 o ) + 4] , 

h 22 = 2G_ csc^^vi [_ 1? CQS + 12 C o S 20 - 3 COS 30 ( 30 ) 

—2 cos(0 — 20 o ) — 3 cos(30 — 20 o ) + 12 cos 20 o 

+6 cos 2(0o - 0) + 6 cos 2(0 + O ) - 15cos(0 + 20 o ) + 4] , 

h 12 = h 21 = jj§r3pv% esc 2 0o sin 2 0/2 [2 sin - sin 20 
- sin 20o + sin 2(0o - 0) + 2 sin(0 + 20 o )] , 

so that the expected strain amplitude h ~ (h 2 x + h\ 2 + 2h\ 2 ) 1 / 2 turns out to be 

2G 

h = p-j/^o csc2 0o {2 [59 cos 2(0o - 0) - cos 0(54 cos (20 o ) + 101)] cos 2 O 
-9 cos (30 - 40 o ) - 9 cos (30 - 20 o ) + 95 cos 20 o + 9 cos40 o - 

sin0[lOlsin20o + 27sin40 o ] + lO6} 1/2 , (31) 

which, as before, strictly depends on the initial conditions of the stellar encounter. A remark is in order at this point. 
A monochromatic gravitational wave has, at most, two independent degrees of freedom while eq. (|30[) seems to show 
three independent parameters associated to the amplitude of the corresponding wave. This is not true since, as usual 
in the TT gauge, we have h + and h x being h + = hu + h 22 and h x = hu + h>2i ■ For details, see 0]. 

As an example, the amplitude of gravitational wave is sketched in Fig. 2 for a stellar encounter close to the Galactic 
Center. The adopted initial parameters are typical of a close hyperbolic impact and are assumed to be b = 1 AU and 
t>o = 200 Kms -1 , respectively. Here, we have fixed Mi — M 2 = 1.4M . 

In the following, we give an estimate of the rate of stellar encounters on hyperbolic orbits in some interesting 
astrophysical conditions as a typical globular cluster or towards the bulge of our Galaxy. 

Let us consider a generic star cluster with mass density profile p(r) composed, for the sake of simplicity, by 
stars of equal mass m moving with relative velocities v Te \(r). Then, the event rate T, that is the number of interacting 
stars on hyperbolic orbits per unit time, is given by 



r = 



f 


~p(r)~ 


/Ro 


m 



2 

,2 



a(r)v re i(r)Airr dr . (32) 



Here, p(r) is the cluster density profile (which we assume to follow a Plummer model Q), R its radius and cr(r ) ~ 
7r[6i(r) 2 — 62 (r) 2 ] the typical cross section of the hyperbolic encounter at distance r from the cluster center (see [ll| for 
details). In evaluating the rate, we are considering only those hyperbolic encounters producing gravitational waves, 
for example, in the LISA range, i.e. between 10 -4 and 10~ 2 Hz (see e.g. [HI)- The integral on the right hand of the 
previous equation can be approximately solved and put in the form 



5.5 x 10~ 10 



10 5 M„ / V lOkms 



For a typical globular cluster (GC) around the Galactic Center, the expected event rate is of the order of 2 x 10 9 
yrs -1 which may be increased at least by a factor ~ 100 if one considers the number of GCs in the Galaxy. If the 
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Figure 2: The gravitational wave-form is sketched as a function of the polar angle (j> for some values of both the impact 
parameter and velocity. In particular, we have fixed Mi = M2 = 1.4Mq. M2 is considered at rest while Mi is moving with 
initial velocity vo = 200 Kms -1 and an impact parameter 6=1 AU. The distance of the source has been assumed to be r — 8 
kpc. As expected, the wave-form has a maximum in correspondence of the peri-astron distance. 



stellar cluster at the galactic center is taken into account and assuming M ~ 3 x 10 6 M©, v ~ 150 km s _1 and R ~ 
10 pc, one expects to have ~ 10~ 5 open orbit encounters per year. If a cluster with total mass M ~ 10 6 M Q , v ~ 150 
km s _1 and R ~ 0.1 pc is considered, an event rate number of the order of unity per year is obtained. These values 
could be realistically achieved by data coming from the forthcoming space interferometer LISA. 

In this letter, the gravitational wave emission on hyperbolic stellar encounters has been analyzed. In particular, 
we have taken into account the expected luminosity and the strain amplitude of gravitational radiation produced in 
a tight hyperbolic impact where two massive objects of 1.4M© closely interact at an impact distance of 1AU. Due 
to the high probability of such encounters inside rich stellar fields (e.g. globular clusters, bulges of galaxies and so 
on), the presented approach could highly contribute to enlarge the classes of gravitational wave sources (in particular, 
of dynamical phenomena capable of producing gravitational waves). In particular, a detailed theory of stellar orbits 
could improve the statistics of possible gravitational wave sources. 
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